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Cosmological scaling solutions, which give rise to a scalar-field density proportional to a back- 
ground fluid density during radiation and matter eras, are attractive to alleviate the energy 
scale problem of dark energy. In the presence of multiple scalar fields the scaling solution can 
exit to the epoch of cosmic acceleration through the so-called assisted inflation mechanism. We 
study cosmological dynamics of a multi-field system in details with a general Lagrangian density 
p = J^"^-^ Xi(/(Xie'^''^'), where Xi — — (V0i)^/2 is the kinetic energy of the i-th field ^i, Ai is a 
constant, and g is an arbitrary function in terms of Yi = XiC'^^'^'- . This covers most of the scalar- 
field models of dark energy proposed in literature that possess scaling solutions. Using the bound 
coming from Big-Bang-Nucleosynthesis and the condition under which the each field cannot drive 
inflation as a single component of the universe, we find the following features: (i) a transient or 
eternal cosmic acceleration can be realized after the scaling matter era, (ii) a "thawing" property 
of assisting scalar fields is crucial to determine the evolution of the field equation of state lu^, and 
(iii) the field equation of state today can be consistent with the observational bound < —0.8 in 
the presence of multiple scalar fields. 



I. INTRODUCTION 

The constantly accumulating observational data con- 
tinue to confirm the existence of dark energy responsible 
for cosmic acceleration today [l|. The cosmological con- 
stant, whose equation of state is w = —1, has been fa- 
vored by the combined data analysis of supcrnovae la , 
cosmic microwave background |3|, and baryon acoustic 
oscillations Meanwhile, if the cosmological constant 
originates from a vacuum energy associated with parti- 
cle physics, its energy scale is enormously larger than 
the observed value of dark energy (pde ~ 10"'^'' GeV''). 
Hence it is important to pursue an alternative possibility 
to construct dark energy models consistent with particle 
physics. 

Scalar-field models such as quintessence 0, @ and k- 
essence have been proposed to alleviate the above 
mentioned problem. In general the energy density of a 
scalar field (f> dynamically changes in time, so that its 
value around the beginning of the radiation-dominated 
epoch can be much larger than the dark energy density 
today. One of such models is quintessence with an expo- 
nential potential V{(j)) = Vbe"^'"'^ where A is a con- 
stant and K — VSttG with G being gravitational constant 
(see Ref. for the classification of cosmological dynam- 
ics and also Refs. [llj for early related papers). In fact, in 
higher-dimensional gravitational theories such as super- 
string and Kaluza-Klcin theories, exponential potentials 
often appear from the curvature of internal spaces asso- 
ciated with the geometry of extra dimensions (so called 
"modulus" fields) [l^. Moreover it is known that expo- 
nential potentials can arise in gaugino condensation as 
a non-perturbative effect [l3l | and in the presence of su- 
jJCT gravity corrections to global supcrsymmetric theories 

The quintessence with an exponential potential 
V{4>) — Voe~^'^'^ gives rise to two distinct fixed points in 



the flat Friedmann-Lemaitre- Robertson- Walker (FLRW) 
background [§]: (a) the scaling solution, and (b) the 
scalar-field dominated solution. If the slope A of the po- 
tential satisfies the condition A > ^3(1 + w/), where Wf 
is the equation of state of a background fluid, then the 
solutions approach the scaling attractor characterized by 
a field density parameter fi^ = 3(1 -I- Wf)/\^ . Even if 
the field energy density is initially comparable to the 
background fiuid density p/, the field eventually enters 
the scaling regime in which is proportional to pf. This 
is attractive to alleviate the fine-tuning problem of the 
energy scale of dark energy. However the scaling solu- 
tion needs to exit from the matter era to the epoch of 
a late-time cosmic acceleration. The scalar-field domi- 
nated solution (ri0 = 1) can be an accelerated attractor 
for A < \/2, but this is incompatible with the condition 
A > ^3(1 -|- Wf) required for the existence of scaling so- 
lutions. Hence the scaling solution cannot be followed by 
the scalar-field dominated solution responsible for dark 
energy. 

There are a number of ways to allow a transition from 
the scaling regime to the epoch of cosmic acceleration. 
One of them is to introduce a single-field potential that 
becomes shallow at late times, e.g., V{(j)) = cie~'*""^ + 
c^e-^""'' with A > y/3{l + Wf) and ^ < ^2 [H] (s§e 
Ref. [l^ for the classification of dynamics and Refs. 
for related works). For this double exponential potential 
the field equation of state lu^ of the final attractor is given 
by = — 1 -I- p^/3. In order to satisfy the observational 
constraint if^ < —0.8 (isj today, we require that p is 
smaller than the order of 1. If the exponential potential 
originates from particle physics models then the slope p is 
typically larger than 1, which is difficult to be compatible 
with the condition for cosmic acceleration. 

Another way is to consider multiple scalar fields with 
exponential potentials, e.g., V{(j>u4'2) = r-.p^^^'^'t'^ 
C2e-^2'='?^2 (see also Refs. UM)- 
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In fact such po- 
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tentials arise from the compactification of higher dimen- 
sional theories to 4-dimensional space-time. It is known 
that the phenomenon called assisted inflation [l^l occurs 
for the multi-field exponential potential, even if the in- 
dividual field has too steep a p otential to lead to cosmic 
acceleration (see also Refs. [23[)- For the sum of steep po- 
tentials satisfying the condition A, > a/2 (i = 1, 2, • • • , n), 
the multiple fields evolve to give dynamics matching a 

single-field model with Aeft = (ELi V^-)"^^^ < ^2 
[2^ . Since the conditions Xi > are mostly satisfied 
for the models motivated by particle physics, this coop- 
erative accelerated expansion is attractive for both infla- 
tion and dark energy. If we apply this scenario to dark 
energy, the scaling radiation and matter eras can be fol- 
lowed by the epoch of assisted acceleration as more fields 
join the scalar-field dominated attractor with an effective 
equation of state = —1 + [20 |. 

The scaling solution arises not only for quintessence 
with an exponential potential but also for more general 
scalar- field models with the Lagrangian density X), 
where X — —g^^ d^(j)du(t> = — (V(?!))^/2 is a kinetic term 
of the field 0. Here g^^ is a metric tensor with the nota- 
tion (—,-!-, -f, -f ). It was found in Refs. i25i] that the 
existence of scaling solutions restricts the form of the La- 
grangian density to be = Xg(Xe^'^^, where A is 
a constant and g is an arbitrary function in terms of y = 
Xe^'^ (here we use the unit ~ 1). The quintessence 
with an exponential potential ~ X — ce~^'^) corre- 
sponds to the choice g = 1 — c/Y, whereas the choice 
g = —1 + cY gives rise to the dilatonic ghost condensate 
model: p = —X + ce^'^X^ [2^ (which corresponds to the 
string-theory motivated generalization of the ghost con- 
densate model proposed in Ref. [1^). The tachyon La- 
grangian density p = —V{(j))^/1 — 2X with V{(l)) oc cf)'^ 
[27| also follows from the above scaling Lagrangian by a 
suitable field redefinition [2^. 

For the multi-field scaling Lagrangian density p = 
X;r=i^45(^»e-^"^') it was shown in Ref. ^ that as- 
sisted inflation occurs with the effective slope Aeff = 

irrespective of the form of g. Hence 
one can expect that the scaling solution is followed by 
the assisted acceleration phase for such a general La- 
grangian. If we consider loop or higher-order derivative 
corrections to the tree-level action motivated from string 
theory (such as e^'^X'^), the constant A is typically of the 
order of unity [sO]. In the single-fleld case this is not 
compatible with the condition for cosmic acceleration. It 
is of interest to see how the presence of multiple fields 
changes this situation. 

In this paper we shall study cosmological dynamics 
of multiple scalar fields with the Lagrangian density 
p = J27=i -^i 9i-^i^^''^')- We are interested in the case 
where the scaling radiation and matter eras induced by 
a field 01 are followed by the dark energy dominated 
epoch assisted by other scalar fields. For the two-field 
quintessence with exponential potentials a similar analy- 
sis was partially done in Ref. [20j , but we shall carry out 



detailed analysis by taking into account bounds coming 
from Big-Bang-Nucleosynthesis (BBN) and supernovae 
observations. In particular the evolution of the field equa- 
tion of state will be clarified in the presence of two and 
more than two fields. We also investigate cosmological 
dynamics for the multi-field dilatonic ghost condensate 
model as an example of k-essence models. 

This paper is organized as follows. In Sec.[TT]we present 
the dynamical equations for our general multi-field La- 
grangian density without specifying any form of g. In 
Sec, mil we derive the fixed points that correspond to the 
scaling radiation/matter solutions and the assisted field- 
dominated attractor. In Sees. IIVI and |V] we study the 
multi-field cosmological dynamics for quintessence with 
exponential potentials and the dilatonic ghost condensate 
model, respectively. Sec. IVIl is devoted to conclusions. 



II. DYNAMICAL SYSTEM 

Let us first briefly review single-field scaling models 
with the Lagrangian density p((j),X). The existence of 
cosmological scaling solutions demands that the field en- 
ergy density = 2Xp,x ~ P, where p^x = dp/dX, is 
proportional to the background fiuid density pf. Under 
this condition the Lagrangian density is restricted to take 
the following form in the flat FLRW background ^ [H] 

p{cl,,X)^Xg{Xe^^), (1) 

where A is a constant and g is an arbitrary function in 
terms of y = Xe^"^. The La grangian density ([T|) is valid 
even in the presence of a constant coupling Q between 
the field <f> and non-relativistic matter and also in the 
presence of a Gauss-Bonnet (GB) coupling between the 
field and the GB term^ [3l|. In the following we do not 
take into account such couplings. Throughout this paper 
we use the unit = SttG = 1. 

The field density parameter for scaling solutions is 
given by fi^ = 3(1 -I- Wf)p^x/^ [1^, HI], where Wf is 
the fluid equation of state. If the field enters the scal- 
ing regime during the radiation era, the BBN places the 
bound il^ < 0.045 at the 2a confidence level |33|. This 
then gives the constraint \^ /p^x > 88.9. 

Besides scaling solutions, there is a scalar-field domi- 
nated point (f20 = 1) with the equation of state = 
-1 + AV(3p,x) m, This can be used for dark en- 
ergy provided that < —1/3, i.e. X^/p,x < 2. Un- 
fortunately this condition is incompatible with the con- 
straint coming from the BBN. Hence the scaling solution 
does not exit to the scalar-field dominated solution in the 
single-field scenario. 



It is also possible to obtain a generalized form of the scaling 
Lagrangian density even when the coupling Q between ij> and 
non-relativistic matter is field-dependent |32|| . 
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If we consider multiple scalar fields <j>i (i = 1, 2, • • • , n) 
with the Lagrangian density 



(2) 



the scaling solution can be followed by the accelerated 
scalar-field dominated point through the assisted infla- 
tion mechanism. Even if the individual field does not sat- 
isfy the condition for inflation, the multiple fields evolve 
cooperativelyto give dynamics matching a single-field 
model with 1291 



f 

^eff 



n 



i=l 



A? 



(3) 



Since Acff is reduced compared to the individual A^, this 
allows a possibility to exit from the scaling matter era to 
the regime of cosmic acceleration. 

In addition to the n scalar fields with the Lagrangian 
density ([2]) we take into account radiation (energy den- 
sity Pr) and non-relativistic matter (energy density Pm)- 
In the flat FLRW space-time with a scale factor a they 
obey the usual continuity equations pr + 4,Hpr = and 
Pm + iHpm = 0, respectively, where a dot represents a 
derivative with respect to cosmic time t and H = a/a is 
the Hubble parameter. The pressure p^- and the energy 
density p^. for the j-th scalar field are given, respectively, 
by 

p^. = X,g{Y,), (4) 
p^. = 2X,p,x. - - X, [g{Yi) + 2Y,g'{Y,)] , (5) 

where a prime represents a derivative with respect to Yt . 
These satisfy the continuity equation 



3H( 



+ =0, 



which corresponds to 



cj), + SHA(Y,)p^xA^ 

+KX^ {1 - A{Y,MY,) + 2Y,g'{Y,)]} = , 

where 

A{Y,) EE [g{Yi) + 5y,g'(y,) + 2Y^g"{Y,)]-' . 
The Fricdmann equations are 



H 



E 



Pm J 



(6) 

(7) 
(8) 

(9) 
(10) 



i=l 



In order to derive autonomous equations we define the 
following quantities 



\/Pr 



(11) 



where the quantity defined in Eq. ([2]) can be expressed 
as 

Y^ = xj/yf . (12) 
We also introduce the field density parameters 



[g{Y,) + 2Y,g'{Y,)] , n^ = J2^^,- 



From Eqs. ^ and ^ it follows that 

Pm 



nr. 



(13) 

(14) 
(15) 



where 0,^ = is the density parameter of radiation. 

Using Eqs. ([6]) and p5|) . we obtain the autonomous 
equations 



dx^ Xi 

div y 

V6 



3 -t- 3 ^ x\g{Yi) ^u^ - VeAiX, 



+ ^A(y,) [a,% - V6{5(K0 + Y,g\Y,)], 



AN ~ 2 

Au u 
A/V 2 



3-f-3^x2<7(y,)+u2- V6A, 

4=1 

n 



(16) 
(17) 

(18) 



where = In (a) . The field equation of state w^p^ of the 
i-th field, the total field equation of state w^, and the 
effective equation of state Wes of the system are given, 
respectively, by 



p^, g{Y,) + 2Y,g'{Yi) 



En 
2=1 - 



Er=iP0> Er=i^i[5m)+2i^.3'(^4)] 



(19) 



,(20) 



(21) 



III. FIXED POINTS OF THE SYSTEM 

Let us derive fixed points for the autonomous equa- 
tions (fT6|) - p8)) . In particular we are interested in the 
scaling solution and the scalar-field dominated solution. 
For these solutions the variables yt do not vanish. Set- 
ting du/dN = in Eq. HH]), it follows that = 
1 — 3 ^"^]^ g(yi) or M = 0. The former corresponds 
to the solution in the presence of radiation, whereas the 
latter to the solution without radiation. In the following 
we shall discuss these cases separately. 
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A. Radiation-dominated scaling solution 

Plugging = 1 - 3J2i=i^^9iYi) into Eqs. p6)) and 
(fT7|) . the fixed point for the i-th field (with 7^ and 
A(Kj) 7^ 0) satisfies 

A,a;, = ^ = ^/6[g(K0 + y,g^(y,)] ^ ^^2) 



ff(r,) + 2r,5'(y,) 



which gives 



(23) 



From Eq. p9|) the field equation of state for the j-th field 
is 



1/3, 



(24) 



which means that p^. is proportional to p,-. Using 
Eqs. (Uni) and together with p^x, ^ 9(^1) +yig' {Yi), 
we have 



4p.x. 
A? 



(25) 



If all n scalar fields are in the scaling regime, Yi are 
the same for all i {Yi = Y) from Eq. ([221) and hence 
p^Xi (* ~ Ij 2, • • • , 71) take a common value p_x = g{Y) + 
y(7'(F) with an effective single-field Lagrangian density 
p = Xg{Y). Then the total field density is given by 



4p,x 



(26) 



where AoS is defined in Eq. We are interested in 

the case where one of the fields, say is in the scaling 
regime in the deep radiation era, while the energy densi- 
ties of other fields are suppressed relative to that of 0i. 
In the BBN epoch we have the following constraint from 
Eq. jlll): 



4p, 



Xi 



A? 



< 0.045 



A. 



> 88.9. 



(27) 



For a given model, i.e. for a given form of g, the vari- 
ables xi and yi are determined by solving Eq. ((22|) . If 
the scalar fields with i 7^ 1 join the scaling solution at 
the late epoch of the radiation era, the total field density 
fl^ tends to increase according to Eq. ([26)) with the de- 
crease of Aeff. If the slope A2 of the second scalar field 
(j)2 that joins the scaling solution is of the order of 1, the 
field density (|26p can be as large as ft^ = 0.1-1. It is not 
preferable for many fields with low Xi to join the scaling 
solution during the radiation era in order to avoid that 
ilcj, exceeds 1. This can be avoided if the field densities 
fl^. (i ^ 1) are much smaller than the radiation density. 



Matter-dominated scaling solution and assisted 
scalar-field dominated point 



the following equations 

n 
i=\ 

KQ^, = V6[giY,) + Y,g'{Y,)]x, . 
From Eqs. ([13|), ([HI), ([28|) and ([29|) it follows that 



, =f]x2.g(yo = -i + ^A.x 



(28) 
(29) 

(30) 



i=l 



Since w^.Vl^^ ~ xfg{Yi) we have 



(31) 



In the case of a single field (f>i , this equation gives w^- ~ 
or rj^. ~ 1. The former corresponds to the scaling solu- 
tion along which p^. is proportional to the matter density 
Pm, whereas the latter is the scalar- field dominated solu- 
tion. 

If all n scalar fields are on the fixed points characterized 
by the condition (|30p , it follows that w^-^ = . . . = 7^;^ . — 
■ ■ ■ ~ w^^ = W(f, and hence Yi = ■ ■ ■ = Yi — ■ ■ ■ ~ Yn = Y 
from Eq. ([TO]) . In this case one has either = or 
fl^ = 1 from Eq. ([31]). Equation ([30|), which holds for 
the each scalar field, reduces to the single-field system 

/fi 

w^^x^g{Y) = -l + ^X,sx, (32) 

where x — XiXi/X^s- The effective single-field La- 
grangian density is given hy p = XgiY) with fi^ ~ 
x^ [g{Y) + 2Yg'\Y)]. We also note that Eqs. ([28|) and 
(|29|l reduce to the following effective single-field forms: 



3-K3a;2.g(r) = \/6Aoff x , 
Aoffri^ = \/Qp^xx , 



(33) 
(34) 



where p^x = 9{Y) +Y g' (Y). 

In the following we shall discuss the matter-dominated 
scaling solution and the assisted field-dominated solution, 
separately. 



1. Matter- dominated scaling solution 

If the i-th scalar field is in the scaling regime during 
the matter-dominated epoch, i.e. w^^ = 0, it follows from 
Eqs. ([HI) and ^ that A^x, = ^6/2 and 



From Eq. ([T5)) we obtain 



(35) 



In the absence of radiation (u ~ 0) the fixed points for 
the i-th field corresponding to j/i 7^ and A{Yi) ^ obey 



(36) 
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More generally the field density parameter in the pres- 
ence of a perfect fluid with an equation of state Wf is 
given by f}^, = 3(1 + Wf)p,xJX^, 

If all n scalar fields are in the scaling regime, then they 
can be described by an effective single-field system with 
= and 



X 



^cff 



(37) 



This scaling solution is stable for A^jf > 3p^x [2£ 



2. Assisted field- dominated point 

Besides the matter scaling solution discussed above, 
there is another fixed point that can be responsible for 
the late-time acceleration. In the single-field case the 
solutions do not exit to the accelerated field-dominated 
point (fJ^. = 1) from the scaling matter era, because 
the scaling solution is stable for il^. = 3p^Xi/Af < 1. 
However the presence of multiple scalar fields allows this 
transition. 

Since il^ = 1 in Eq. ((34|) for the scalar-field dominated 
point with n multiple fields, it follows from Eq. ([5^ that 



1 



3p. 



X 



(38) 



This fixed point can be responsible for the late-time ac- 
celeration (ui0 < —1/3) for Apff < 2p^x- Moreover it 
is stable under the condition < 3p^x [2§| (which is 
opposite to the stability of the scaling matter solution) . 
Using the relations w^fi^. = xjgiY) and = x^g{Y), 
we find 



(39) 



We shall study the case in which one of the fields has 
a large slope Ai (» 1) to satisfy the BBN bound ((27l) 
and other fields with A^ — 0(1) join the scalar-field dom- 
inated attractor at late times. Since the joining of 
such multiple scalar fields reduces Aoff it should be possi- 
ble to give rise to sufficient cosmic acceleration through 
the assisted infiation mechanism, even if the individual 
field cannot be responsible for the acceleration. 

For a given model one can derive Yi (for the field 4>i) 
that corresponds to the scaling solution during radiation 
and matter eras by solving Eqs. (^5)1 and ([55]) . respec- 
tively. The field density parameters fi^j in these epochs 
are given by Eqs. (^5]) and respectively. The as- 

sisted field-dominated solution corresponds to 



6[g{Y)+Yg'(Y)Y 
g{Y) + 2Yg'{Y) 



A; 



(40) 



which comes from by combining Eqs. ([55]) and ([M]) with 
= 1. By solving this equation for a given form of 



g{Y), we obtain the field equation of state ((38|) and also 
X = Aeff/(V6p,x) from Eq. 

In subsequent sections we shall consider two models: 

(i) quintessence with multiple exponential potentials, and 

(ii) the multi-field dilatonic ghost condensate model (one 
of k-essence models). In our numerical simulations we 
identify the present epoch (the redshift z = 0) to be 
il0 = 0.72 with the radiation density in the region 7.0 x 
10-^ <Vlr< 1.0 X 10~"*. 



IV. QUINTESSENCE WITH MULTIPLE 
EXPONENTIAL POTENTIALS 



The single-field quintessence with an exponential po- 
tential corresponds to the Lagrangian density p = X — 
i.e. the choice g{Y) = 1 — c/y in Eq. In 



-X<p 



the following we shall consider the Lagrangian density 
([2]) of n scalar fields with the choice g{Yi) = 1 — Ci/Yi 
(i = l,2,---,n). 

Since p^Xi = 9{Yi) + Yig'{Yi) = 1 in this model the 
scaling field density fl^. during the radiation and matter 
eras is given by flif,. — A/Xf and Q^. = 3/A?, respectively 
[see Eqs. (pSj) and ((551) ]. Below we discuss the case in 
which one of the scalar fields, (jfi, is in the scaling regime 
during most of the radiation and matter eras and other 
fields eventually join the assisted scalar-field dominated 
attractor with given by Eq. ((55)) . Then the BBN 
bound ((?7|) gives 



Ai > 9.42. 



(41) 



Under this condition, the scaling field density il^^ dur- 
ing the matter-dominated epoch is constrained to be 
il^j < 0.034. If other fields join the scaling regime in the 
radiation (matter) era, the field density increases from 

n^, = 4/ XI (n^, = 3/ A?) to = 4/xlg {n^ = s/a^^,). 

This is possible provided that the slopes of the joining 
scalar fields satisfy the conditions A; ^ 1. Meanwhile, 
if Ai are of the order of 1, this leads to a large density 
parameter fl^ that is comparable to unity. In what fol- 
lows wc focus on the case in which the fields with slopes 
Ai — 0(1) {i > 2) enter the regime of the assisted cosmic 
acceleration preceded by scaling solutions induced by (f)i . 
It is convenient to introduce the following variable 



From Eq, 
for the field 



y, = ^y,. (42) 

the radiation-dominated scaling solution 
corresponds to 



'2V6 2^/t 
3A7' 3A7 



Yi = 2ci 



(43) 



This is followed by the matter-dominated scaling solu- 
tion, satisfying 



V6 j/6 
2Ai' 2Ai 



= ci . (44) 
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The assisted field-dominated point corresponds to the 
single-field potential = ce"^"""^, i.e. g{Y) = 1-c/Y 
with Y = I'lp-. From Eqs. (gOl) and ([34|) this is charac- 
terized by the fixed point (where we define y = y/cy): 



1 - 



6 



Y 



(45) 



For the i-th field we have that Xi ~ (Xcff/Xi)x and yi = 
x,/VY. 



A. Two fields 

First let us consider the case of two scalar fields <j>i and 

In Fig. [T] we plot the evolution of the background fluid 
density pf = Pr + Pm and the field densities p^^ , ver- 
sus the redshift z = ao/a — 1 (oq is the present value of 
a) for Ai = 10 and X2 — 1.5. We choose three different 
initial conditions for 0i. The case (i) corresponds to the 
exact scaling solution starting from the fixed point (j43p . 
along which Vl^^ = = 0.04 and f^^, = S/Af = 0.03 
during radiation and matter eras, respectively. Finally 
the system enters the epoch in which the energy density 
of the second field (^1 dominates the dynamics. Fig- 
ure [T] shows that the field 0i eventually joins the scaling 
regime both for the initial conditions (ii) p^^ ~ Pf and 
(iii) Pf- Thus the cosmological trajectories con- 

verge to a common scaling solution for a wide range of 
initial conditions. 

In Fig. [1] we find that the second field density is 
almost frozen after the initial transient period. In order 
to understand this behavior we introduce the ratio r.i 
between the kinetic energy 0f /2 and the potential energy 
V,{<j)i) = c.e"^''^- of the i-th field: 



2V^ 



(46) 



which is related to the quantity Yi via = Yi/ci. Taking 
the derivative of r.i with respect to N , it follows that [3^ 



d In ri 



6[A,(t)-l] 



A,(t) = A, 



3(l + ^«0j' 



For the scaling field (pi one has w^-^ = ^/ and fJ^^ = 
?>{\ + Wf)/X\ 0, where is the equation of state of 
the background fluid. This means that Ai(t) = 1, so 
that the ratio ri = Yi/ci remains constant. In fact, from 
Eqs. (|43|) and (|44|) . one has ri = 2 and ri ~ 1 during the 
radiation and matter eras, respectively. This reflects the 
fact that the scaling field has a kinetic energy with the 
same order as its potential energy. 

The field ^2 joining the assisted attractor at late-times 
satisfies A2 = ^ Ai and Vl^^ <C fi^i at the early 

stage of the radiation era, so that ^2{t) ^ 1 initially 
(unless is unnaturally close to —1). At this stage 



a 
o 



0) 






12 10 8 6 4 
logio(1+z) 



2 -2 



Figure 1: Evolution of the energy densities p/ = pr + pm, P0i, 
and p<j42 versus the redshift z for Ai = 10 and A2 = 1.5 in the 
two-field quintessence with exponential potentials. We choose 
three diflterent initial conditions for (j>i: (i) xi = 2\/6/(3Ai), 
yi = 2a/3/(3Ai), (ii) xi = 0.99, yi = 0.12, and (iii) xi = 
1 X 10^^, yi = Ix 10~^, while other initial conditions are fixed 
to be 2:2 = 1 X 10"^\ ^2 = 2 X 10~^^, and = 4 x 10"^°. 
In the case (i) the field is in the scaling regime from the 
beginning. 



the ratio r2 decreases rapidly as cx e~^^ according to 
Eq. (j47|) , see Fig. [2l In the region r2 < 1 the field 02 is 
almost frozen with nearly constant p^g • As r2 decreases, 
grows and approaches —1. This leads to the 
growth of A2 (t) . As we see in Fig. [2] the ratio r2 starts 
to increase after A2{t) becomes larger than 1. When r2 
grows to the order of 1, the field 02 begins to evolve to 
join the assisted attractor given by Eq. (gS]). 

The mass squared for the i-th scalar field is given by 
mf = d2y,(4)/d02 = A2y,(0,). The energy density of 



3Hl 



starts to dominate around the present epoch, so that 
^2(02"'') (t^*^ subscript "0" represents present val- 



(47) ues). Then the mass of 02 can be estimated as 



m2(0f ) 



X2H0 



(48) 



Recall that A2 needs to be of the order of 1 to realize a 
stable assisted attractor satisfying the condition A^^j < 
3. Hence the mass m2 is as small as today. In the 
numerical simulations of Figs. [T] and [5] the field 02 is 
almost frozen with the mass P5|) for the redshift 1 < z < 
10^^ (during which the condition r2 <C 1 is fulfilled). 

Even if the field 02 is rapidly rolling at the initial stage 
of the radiation era such that r2 ^ 1, it enters the regime 
in which 02 is nearly frozen (w^j ~ —1) prior to the 
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Figure 2: Evolution of r2 and A2 — 1 for Ai = 10 and A2 = 1.5 
in the two-field quintessence with exponential potentials. The 
same initial conditions are chosen as in the case (i) of Fig. [1] 
Initially the quantity r2 decreases as r2 oc e~®^ because A2 ~ 
0. The ratio r2 starts to increase after A2 becomes larger than 
1. 



matter-dominated epoch. In Fig.[3]the evolution of p^^ is 
plotted for three different initial conditions of X2 with y2 
fixed. The dominance of the field kinetic energy relative 
to its potential energy corresponds to r2 ^ 1 and w^^ ~ 
1, which results in the rapid decrease of il^j to reach 
the regime ^2{t) <^ 1. Even if r2 ^ 1 initially, the 
decrease of r2 in the regime A2(i) ^ 1 is so fast (oc e^^^) 
that the field (j)2 eventually enters the frozen regime with 
~ —1. For the initial conditions satisfying r2 ^ 1 
the field 02 is almost frozen from the beginning, so that 
/902 is nearly constant until recently. 

If we change the initial conditions of j/2 associated with 
the field potential, this leads to the modification of the 
epoch at which the field (j)2 dominates at late times. This 
comes from the fact that the density p^^ during which (j)2 
is nearly frozen is sensitive to the choice of its initial po- 
tential energy. Thus the evolution of the field ^2 depends 
on its initial potential energy but not on its initial kinetic 
energy. 

Figure [4] illustrates the variation of w^, w^^, w^^, and 
WcS for Ai = 10 and A2 = 1.5 with the same initial 
condition as in the case (i) of Fig. [TJ The equations of 
state and w^-^ are similar to the effective equation 
of state Wcff during radiation and matter eras, but the 
deviation appears at low rcdshifts. The field (/)2 is almost 
frozen around w^^ ~ ~1 after the initial transient period, 
but it begins to evolve for z < 0{1). 



Figure 3: Evolution of p/ , p^j, and p^^ for Ai = 10 and A2 = 
1.5 in the two-field quintessence with exponential potentials 
with three different initial conditions for the kinetic energy of 
02: (i) a-2 = 1 X 10"2\ (u) X2 = lx 10~", and (iii) X2 = 0.05. 
Other initial conditions are chosen to be j/2 = 2 x 10^^*, 
XI = 2^/6/(3Al), yi = 2^/(3Ai), and n„ = 4 x 10"^^ 
The field 02 enters the regime with a nearly constant 
independent of its initial kinetic energy. 

From the definition of in Eq. ([^0]) we have 

if^ = 7^ (w^jfi^i -I- Wtji^n^,^) . (49) 

Note that w^-^ w and w^^ ~ —1 around the end of 
the matter-dominated epoch. After gets larger than 
il^j^ , Wrj, begins to be mainly determined by the field 02 , 
i.e. W0 « W02r2^2/^0- As we see in Fig. 2] takes 
a minimum before reaching the present epoch (z = 0), 
which is followed by its increase toward the attractor 
value — —1 + Agg/3. For the model parameters used 
in the numerical simulation of Fig.|4]we have AeS = 1.483, 
which gives = —0.267 at the scalar-field dominated 
attractor. This corresponds to the decelerated expansion 
of the universe. Meanwhile one has w^{z = 0) = —0.62 
and Wcs{z = 0) = —0.45, which means that the transient 
acceleration occurs at the present epoch. Interestingly, 
even without the assisted accelerated attractor, such a 
temporal acceleration can be realized by the presence of 
the thawing field 02. 

Under the BBN bound (gT]) and the condition A2 > ^ 
(i.e. the field 02 cannot be responsible for the accelerated 
expansion as a single component of the universe), the 
equation of state for the late-time assisted attractor 
is not very different from —1/3. Meanwhile the present 
value of Wtj, is smaller than its asymptotic value. For the 
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Figure 4: Evolution of the field equations of state w^, w^^^, 
111^2 1 E'lid the effective equation of state Wcff for Ai = 10 and 
A2 = 1.5 in the two-field quintessence with exponential poten- 
tials. The initial conditions are chosen to be xi = 2\/6/(3Ai), 
yi = 2^/(3Ai), X2 = 1 X 10"^\ y2 = 2 X 10"^-*, and 
= 4 X 10-^°. 
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Figure 5: The field equation of state today versus A2 for 
Ai = 9.43 (solid curve) in the two-field quintessence with ex- 
ponential potentials. If A2 < 1.76 the scalar-field dominated 
point is the final attractor. The condition for cosmic acceler- 
ation (wctii^z = 0) < —1/3) is satisfied even for A2 > \/2. The 
two lines (a) and (b) in the figure correspond to A2 = 1.76 
and A2 — V2, respectively. 



marginal case with Ai = 9.43 and A2 = 1.415 we find 
that ^0(2 = 0) = —0.66 numerically. For increasing A2 
we obtain larger values of w^{z = 0) and Wcs{z — 0), 
as we see in Fig. [5j If we do not impose the condition 
A2 > V2, then 'w^{z = 0) can be smaller than —0.66. 
Note that, when A2 > 1.76 and Ai = 9.43, the scalar-field 
dominated point ceases to be the late-time attractor. We 
have also carried out numerical simulations for different 
values of Ai satisfying the condition Ai > 9.42 and found 
that Wif,{z = 0) and Wcs{z = 0) are insensitive to the 
change of Ai. 



B. More than two fields 

For three scalar fields the total field equation of state 
is given by = {w^^D,^j^+w^2D,^2+w^^n^^)/n4,. If the 
two fields 02 and ^3 with slopes A2, A3 = 0(1) join the 
assisted attractor for z < 0(1), it is possible to obtain 
smaller values of w^{z = 0) and Wcs{z ~ 0) relative to 
the two-field case. 

In Fig. [6] we plot the evolution of w^, w,f,- {i ~ 1,2,3) 
as well as Wcs for Ai = 9.5, A2 = 1.42, and A3 = 2.0. 
The field 0i is in the scaling regime during the radiation 
and matter eras, which is followed by the epoch of cosmic 
acceleration once the energy densities of 02 and 03 are 
dominant. The fields 02 and 03 have been nearly frozen 
(except for the initial transient period) by the time they 
start to evolve for z < 0(1). In the numerical simulation 
of Fig. [6] the energy densities p^^ and are the same 
order when they begin to dominate over the background 
fluid density. In the numerical simulation of Fig. [S] the 
field equation of state today is found to be w^{z = 0) = 
—0.76, which is smaller than the minimum value w^{z = 
0) = —0.66 in the two-field case. This comes from the 
fact that the third field with lu^g close to —1 leads to 
smaller values of w^{z = 0). 

For the marginal model parameters Ai ~ 9.43, A2 = 
A3 — 1.415, which satisfy the conditions Ai > 9.42 and 
A2,A3 > V2, we find that ■w^{z = 0) = —0.83, provided 
the fields 02 and 03 exit from the frozen regime almost 
at the same time. If either 02 or 03 begins to evolve 
much later than another, then 'w^{z = 0) tends to be 
larger. This case is not much different from the two-field 
scenario for estimating the value of w^{z = 0), although 
the scalar-field dominated attractor is different. In the 
three-field scenario the equation of state for the at- 
tractor can be as small as ~ —0.6 for A2, A3 = 0(1) 
so that cosmic acceleration today is not transient. 

In Fig. [7] we show w^{z = 0) versus A^ {i > 2) for 
Ai = 9.43 in the presence of multiple scalar fields. Under 
the condition A^ > the observational bound, Wff,{z = 
0) < —0.8, can be satisfied for three fields. In the case 
of four fields it is possible to satisfy the same bound for 
Aj < 2 (i = 2,3,4). As we add more fields, we obtain 
smaller values of w^{z = 0) and Weg{z = 0). 
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Figure 6: Evolution of w^, w^-^, w^^j '^</>3) ^-nd '''eft for Ai — 
9.5, A2 = 1.42, and A3 — 2.0 in the three-field quintessence 
with exponential potentials. The initial conditions are chosen 
to be xi = 2V6/(3Ai), yi = 2^/(3Ai), X2 = I x 10'^^, 
j/2 = 2 X 10"^^, = 1 X 10-^\ ya = 2x 10"^*, and il^ = 
5 X 10" "\ 



V. MULTI-FIELD DILATONIC GHOST 
CONDENSATE MODEL 

Let us next proceed to the dilatonic ghost condensate 
model with n scalar fields, where the Lagrangian den- 
sity is given by ([2]) with the choice g{Yi) = — 1 -f- aYi 
(^ = l,2,---,n), i.e. P = Etii-^^ + c,e^''*''^?)- The 
coefficients Ci are positive so that the quantum instabil- 
ity problem of the negative kinetic energy {—Xi) can be 
avoided by the presence of the higher-order derivative 
term ae^'-^'Xf 

In this model we have p,Xi = 2Yi — 1 and 



3Y-1 



=a;2(3f, -1), 



where 



(50) 



(51) 



When Yi = 1/2 the equation of state w^. is equivalent to 
— 1. The quantum stabihty of the scalar field is ensured 
for Yi > 1/2 (i.e. > —1), whereas in the region 

Yi < 1/2 the vacuum is unstable against the catastrophic 
particle production of ghost and normal fields [1^ . In the 
following we shall focus on the case % > 1/2. 

From Eqs. ([TB]) and (fTT)) we obtain the following equa- 



Figure 7: The equation of state w,t> today versus Ai {i 7^ 1) 
for Ai = 9.43 in the multi-field quintessence with exponential 
potentials. The slopes Ai for the i-th (i > 2) fields are chosen 
to be the same. In this simulation the fields (f>i (i > 2) enter 
the regime of the assisted field-dominated attractor almost at 
the same time. Under the condition Ai > \/2 (which is shown 
as a thin dotted line in the figure), we have that Wcj, < —0.8 
for more than two fields. 



tions for K^: 

d^ ^ 
dN 



-SYXV&XtXi 



4) + 6 - \/6A,; 



QY,-l 



(52) 



which hold for i ~ 1, 2, • • • , ?i. We will solve Eqs. pi 
(fT8)) and ((52|) in our numerical simulations. 

For this model the solution to Eq. ([^5]) does not exist, 
whereas the solution to Eq. ((35)) is given by = 1 . This 
means that the scaling solution is absent during the radi- 
ation era, while it is present during the matter era. More 
precisely, for the background fluid with an equation of 
state w/, the presence of the scaling solution corresponds 
to the condition (1 - w^)g{Yi) = 2wsY,g'{Yi) [l!]. Solv- 
ing this equation for the present model, we obtain 



1 -3u;/ 



(53) 



For the radiation fluid (wf ~ 1/3) we require that 
Yi 00 for the existence of the scaling solution. If the 
field 01 is in a nearly scaling regime during the radiation 
era, it follows that fl^^ ~ (4/A^)(2Yi - 1). The BBN 
bound fi^, < 0.045 amounts to 



Ai > 9.42 V 2 Yi - 1 



(54) 
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This shows that, under the condition Yi oo, Ai needs 
to be infinitely large. However, as long as we do not 
demand the exact scaling radiation solution, the variable 
Yi can be of the order of unity (as we will see later). In 
such a case the constraint on Ai is not so restrictive. 

The radiation-dominated epoch can be followed by 
the scaling matter era characterized by the fixed point 
{xi,Yi) = (V6/(2Ai),l) with fl^, = 3/Af. The solu- 
tions finally approach the assisted field-dominated point 
satisfying Eq. ([iD|) . i.e. 

where c is the coefficient of the effective single-field La- 
grangian density: p = —1 + cX. In deriving Eq. (j55p we 
have taken the solution with Y > 1/2. From Eq. pS]) 
the field equation of state is given by 

which shows that the late-time cosmic acceleration occurs 
for Aoff < VG/S. The stability of this solution is ensured 
for A2ff < 3(2? - 1), i.e. Kg < ^/3. 

A. Two fields 
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Figure 8: Evolution of S7r, f2m, f2<pi, and Q^^ versus the red- 
shift z for Ai = 40 and A2 = 1 in the two-field dilatonic ghost 
condensate model. The initial conditions are chosen to be 
xi = 2V6/(3Ai), Yi = 5.0, X2 = 1.6 x 10"^", Y2 = 10.0, and 
f7m = 8 X 10"^. AVhile the BBN bound fl^^ < 0.045 is satis- 
fied at z ~ 10^°, the growth of ft^,^ continues by the redshift 
at 2 = 7.2 X 10* with the maximum value fi^j = 0.071. 



We first study cosmological dynamics of the two-field 
ghost condensate model. 

Let us consider the case in which the field 0i initially 
exists around xi ~ 2'\/6/(3Ai) with a finite value of Yi 
satisfying the condition Yi > 1/2 [see Eq. (HH)]. From 
Eq. ((52|) it follows that dYi/dN > 0, as long as xi does 
not depart significantly from 2-\/6/(3Ai). This means 
that the quantity Yi tends to grow during most of the 
radiation era, which also leads to the increase of the den- 
sity parameter fJ^^ = xl{3Yi — 1). This growth of Yi is 
associated with the fact that the radiation scaling solu- 
tion exists only in the limit Yi 00. 

In Fig. [5] we plot one example about the evolution of 
density parameters for Ai = 40 and A2 = 1. This shows 
that ri^j in fact increases during the radiation-dominated 
epoch. In this case we have $7^^ = 0.036 around the 
BBN epoch (z « 10^°), so that the bound n^, < 0.045 
is satisfied. The growth of fi^j ceases around the end of 
the radiation era, because Yi begins to decrease toward 
the scaling matter fixed point at Yi = 1. In order to 
satisfy the BBN bound ([54]) we have numerically found 
that Ai is required to be at least larger than 30. This 
is related to the fact that, even for the initial conditions 
of Yi close to 1/2 around the beginning of the radiation 
era, Yi grows to be larger than 5 at the BBN epoch. 

In Fig. [5] we find that the scaling matter era in which 
rj^j is nearly constant is very short, unlike the multi-field 
quintessence with exponential potentials. This can be 
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Figure 9: Evolution of w^, w^^, w^^j s^nd Woff for Ai = 40 
and A2 = 1 in the two-field dilatonic ghost condensate model. 
The initial conditions are the same as in Fig. |8] 
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Figure 10; The equations of state and w^b today versus A2 
for Ai = 40 in the two-field dilatonic ghost condensate model. 
If A2 < 1.734 the scalar-field dominated point is the final 
attractor. If A2 > the field (j)2 cannot be responsible 

for cosmic acceleration as a single component of the universe. 
The two lines (a) and (b) in the figure correspond to A2 = 
1.734 and A2 = \/6/3, respectively. 



understood as follows. Around the end of the radiation- 
dominated epoch the quantity Yi has already increased 
to a value larger than the order of unity. It takes some 
time for the solutions to reach the scaling matter fixed 
point at Yi = 1. In the numerical simulation of Fig. [8] 
this happens for the redshift at 6.2. Since the solutions 
enter the dark energy dominated epoch for z < 0(1), the 
period of the scaling matter era is short. 

Figure [H] illustrates the evolution of ly^, w^^, w^j, and 
Weft for the same model parameters and initial conditions 
as given in Fig. [51 Initially w^-^ is smaller than Wcs, but 
it grows to the value close to WcH = 1/3 during the radi- 
ation era with the increase of Yi . The evolution of w^-^ 
in Fig. [9] clearly shows that the field 0i does not soon 
enter the scaling matter regime just after the radiation- 
dominated epoch. The field 4>2 approaches the phase 
with Wtjt^ — —1 after the initial transient period. This 
corresponds to Y2 — 1/2 and fi^a ~ 2^2 /2 in Eq- (jSU)) . 
Numerically we find that the late-time cosmological evo- 
lution is practically independent of the initial conditions 
of Y2, but it is sensitive to the initial values of X2 be- 
cause the quantity X2 is associated with the dark energy 
density. 

The epoch at which the field 02 starts to exit from 
the regime w^^ = — 1 depends on the parameter A2. For 
decreasing A2 the redshift Zc at which this "thawing" oc- 
curs gets smaller, which leads to smaller values of and 
Wes today. In Fig. [TO] we plot w^{z = 0) and Wcsiz = 0) 



Figure 11: The equation of state today versus A; (i 7^ 1) for 
Ai = 40 in the multi-field dilatonic ghost condensate model 
(up to five fields with the same slopes Xi, i > 2). In this 
simulation the fields (j>i (i > 2) enter the regime of the assisted 
attractor almost at the same time. Under the condition Ai > 
\/6/3 (which is shown as a thin dotted line in the figure), it 
is difficult to realize ^^(z = 0) < —0.8 even in the presence 
of five scalar fields. 

versus A2 for Ai — 40. In this case the stability of the 
assisted field-dominated point is ensured for A2 < 1.734. 
If A2 > a/6/3 the field 02 cannot drive cosmic accelera- 
tion as a single component of the universe. Under these 
bounds we find that the condition Wcsiz = 0) < —1/3 
for the acceleration today is not satisfied in the two-field 
case. This is intimately associated with the fact that the 
"thawing" of the field 02 occurs quite early (zc > C(10)) 
for A2 > \/6/3, see Fig. [9] This property is different 
from two-field quintessence with exponential potentials 
in which 02 begins to evolve at smaller redshifts for the 
same values of A2. We also note that the present values 
of and Wcff are insensitive to the choice of Ai, as long 
as the condition is satisfied. 

B. More than two fields 

In the presence of more than two scalar fields it is pos- 
sible to obtain smaller values of Wtj, and Wcs today rela- 
tive to the two-field scenario discussed above. In Fig. [11] 
we plot w^{z = 0) versus Ai {i > 2) for Ai = 40. The 
initial conditions for 0i {i > 2) are chosen so that they 
join the assisted attractor almost at the same time. For 
the three-field scenario with A2 = A3 = 0.817 (slightly 
larger than •\/6/3) one has w^{z = 0) = —0.60 and 
Weg{z = 0) = —0.43, so that cosmic acceleration is real- 
ized today. However this case is difficult to be compatible 
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Figure 12: The present equations of state 'w^{z = 0) and 
■ii)off(z = 0) and at the assisted attractor versus the num- 
ber n of scalar fields for Ai = 40 and Ai = 0.817 {i > 2) in 
the multi-field dilatonic ghost condensate model. In this sim- 
ulation the fields <j)i (i > 2) enter the regime of the assisted 
attractor almost at the same time. 



with the observational bound w^{z = 0) < —0.8. Even 
for the five-field case with Ai = 0.817 {i = 2,3,4,5) we 
find that ■^^(z = 0) = —0.70, which is still larger than 
-0.8. 

In Fig.[T2]we plot 'w^{z — 0), lUcff = 0), and at the 
late-time attractor for Ai = 40 and A^ = 0.817 (i > 2). 
This shows that wc require at least 10 scalar fields to 
realize the condition w^{z = 0) < —0.8. Equation ((SS]) 
leads to larger Wcj, at the assisted attractor relative to 
the case of the multi-field quintessence with exponential 
potentials for the same values of Acff. In addition to the 
early "thawing" of assisting scalar fields, this is another 
reason why a large number of fields are required to realize 
small Wcf,{z — 0) close to —1. In Fig. [12] we find that 
w^{z = 0) is almost proportional to u;^ (attractor) for n > 
3. Unless we have many fields such that n > 10, Wtf,{z = 
0) as well as (attractor) are not reduced sufficiently to 
satisfy the observational bound. Of course, if we do not 
demand the condition A^ > \/6/3 {i > 2), it is possible 
to realize w^{z = 0) < —0.8 without introducing many 
fields. 



VI. CONCLUSIONS 

In this paper we have studied cosmological dynamics 
of assisted dark energy for the Lagrangian density ([2]) 
that possesses scaling solutions. This scaling Lagrangian 
density involves many models such as quintessence with 
exponential potentials, dilatonic ghost condensates, and 



tachyon fields with inverse power-law potentials. As long 
as the energy density of a field (with Af » p^Xi ) dom- 
inates over those of other fields, the density parameter 
ri^j remains constant during the radiation and matter 
eras {V,^^ = Ap^Xi/^i and il^^ = 3p^Xi/^i, respectively). 
This property is attractive because the solutions enter the 
scaling regime even if the field energy density is initially 
comparable to the background fluid density. 

In the presence of multiple scalar fields the scaling 
matter era can be followed by the phase of a late-time 
cosmic acceleration as long as more than one field join 
the assisted attractor. The field equation of state for 
the assisted attractor takes an effective single-field value 
ui^ = — 1 + \lff/{3p^x), with Aoff given by Eq. Since 
Aoff is smaller than the slope A^ of the each field, the pres- 
ence of multiple scalar fields can give rise to cosmic ac- 
celeration even if none is able to do so individually. This 
is a nice feature from the viewpoint of particle physics 
because there are in general many scalar fields (dilaton, 
modulus, etc) with the slopes Ai larger than the order of 
unity. 

While the above property of cosmological dynamics is 
generic for the scaling models with the Lagrangian den- 
sity the evolution of as well as $7^ is different 
depending on the forms of the Lagrangian density p. In 
order to see this we have focused on two models: (i) 
canonical fields with exponential potentials, and (ii) mul- 
tiple dilatonic ghost condensates. These correspond to 
representative examples of quintessence and k-essence, 
respectively. 

For the multi-field quintessence with exponential po- 
tentials, the slope Ai for the scaling field is constrained to 
be Ai > 9.42 from the BBN bound. We have numerically 
found that the transient cosmic acceleration today with 
a non-accelerated attractor can be realized after the scal- 
ing matter era. This comes from the thawing property 
of assisting scalar fields that start to evolve only recently 
from a nearly frozen regime characterized by the equa- 
tion of state w,p. ~ —1 {i > 2). Even for the initial 
conditions where the kinetic energies of {i > 2) are 
much larger than their potential energies, we have con- 
firmed that the fields 0, enter the frozen regime by the 
end of the radiation-dominated epoch. In the presence of 
three scalar fields we have found that the total field equa- 
tion of state today can be smaller than —0.8, even if 
each field is unable to be responsible for the accelerated 
expansion as a single component of the universe. 

The multi-field dilatonic ghost condensate model does 
not possess an exact scaling radiation era, although the 
scaling matter era is present. In this model the slope 
Ai of the field 0i is more severely constrained from the 
BBN bound relative to the multi-field quintessence with 
exponential potentials. The fields <j>i [i > 2) enter the 
regime characterized by w^. ~ — 1 during the radiation- 
dominated epoch. However the exit from this regime 
occurs earlier than the multi-field quintessence with ex- 
ponential potentials for the same values of Ai, which gen- 
erally leads to a larger field equation of state today. 
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In the two-field scenario we have found that cosmic accel- 
eration does not occur at the present epoch if the fields 
are unable to give rise to inflation individually. While 
the acceleration today is possible in the presence of more 
than two fields, we require at least 10 fields to satisfy 
the observational bound w^{z = 0) < —0.8 under the 
condition > \/6/3 {i> 2). 

In single-field scaling models with an exit to the late- 
time acceleration (such as the model in Ref. [HI), the 
field equation of state changes from to negative dur- 
ing the transition from the matter era to the accelerated 
epoch. Meanwhile the multi-field dark energy models we 
have discussed in this paper exhibit a rather peculiar be- 
havior of Wfj): it first reaches a minimum and then starts 



to grow toward the assisted attractor (see Figs. |4] and [9]). 
It will be of interest to see whether future high-precision 
observations will detect some signatures for such dynam- 
ics. 
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